Above-room-temperature ferromagnetism in half-metallic Heusler compounds NiCrP, 
NiCrSe, NiCrTe and NiVAs: A first— principles study 
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We study the interatomic exchange interactions and Curie temperatures in half-metallic semi 
Heusler compounds NiCrZ (Z=P, Se, Te) and NiVAs. The study is performed within the framework 
of density functional theory. The calculation of exchange parameters is based on the frozen-magnon 
approach. It is shown that the exchange interactions in NiCrZ vary strongly depending on the Z 
constituent. The Curie temperature, T c , is calculated within the mean field and random phase 
approximations. The difference between two estimations is related to the properties of the exchange 
interactions. The predicted Curie temperatures of all four systems are considerably higher than room 
temperature. The relation between the half-metallicity and the value of the Curie temperature 
is discussed. The combination of a high spin-polarization of charge carriers and a high Curie 
temperature makes these Heusler alloys interesting candidates for spintronics applications. 

PACS numbers: 75.50.Cc, 75.30.Et, 71. 15. Mb 



INTRODUCTION 

An important current problem on the way to the prac- 
tical use of the spin-transport in semiconductor devices 
is the fabrication of the materials that make possible the 
injection of spin-polarized electrons into a semiconductor 
at room temperature 0. One of the classes of systems 
promising to supply materials with necessary combina- 
tion of properties is the half-metallic ferromagnets, i.e., 
ferromagnetic systems where one spin channel is metallic 
and other is semiconducting. The half-metallicity leads 
to the 100% spin-polarization of the electron states at the 
Fermi level Ep 0, El • The half-metallic ferromagnetism 
(HMF) was discovered by de Groot et al. when studying 
the band structure of semi Heusler compound NiMnSb 
. Ishida et al. have proposed that also the full-Heusler 
alloys of the type Co2MnZ, (Z=Si,Ge), are half-metals 
0. Since then a number of further systems were pre- 
dicted to be half-metallic. Among them binary mag- 
netic oxides (Cr02 and FeaO^), colossal magnetoresis- 
tance materials (S^FeMoOg and Lao.ySro.sMnOs) 0, di- 
luted magnetic semiconductors (Gai-^Mn-j-As) and zinc- 
blcndc compounds MnAs and CrAs 

The Heusler alloys form a particularly interesting class 
of materials since they are characterized by a high 
Curie temperature and good crystallographic compatibil- 
ity with the lattices of semiconductors used industrially 
|l0l IllL Il2| . The main body of previous first-principles 
studies was focused on the investigation of the properties 
of the semiconducting gap in the minority-spin channel. 
Galanakis et al. considering NiMnSb have shown that the 
gap arises from the interaction between the ci-orbitals of 
the Ni and Mn atoms leading to the formation of bond- 
ing and antibonding states separated by a hybridizational 
gap 0]. Chioncel et al. have demonstrated by the ex- 
ample of NiMnSb that the electron-magnon interaction 
can lead to the appearance of nonquasiparticle states in 



the half- metallic gap 0] . 

For a number of Heusler alloys it was shown that half- 
metallicity is preserved under tetragonalization of the 
crystal lattice [15| and application of the hydrostatic 
pressure 0]. Mavropoulos et al. studied the influ- 
ence of the spin-orbit coupling on the spin-polarization 
at the Fermi level and found the effect to be very small 
[l?! that corresponds to a small orbital moment calcu- 
lated by Galanakis 

0- 

Orgassa and collaborators have 
shown that the half-metallic gap decreases with increas- 
ing disorder Picozzi and collaborators and Miura 
et al. demonstrated that different types of structural 
defects have different influence on the half-metallicity 
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An important part of theoretical efforts is a first- 
principles design of new half-metallic Heusler alloys. 
Galanakis studied the appearance of half-metallic ferro- 
magnetism in quaternary Heusler alloys |22|. Xing et 
al. have predicted the half metallic ferromagnetism in 
NiCrZ (Z=P, Se, Te) and NiVAs semi Heusler compounds 
00 

Despite a very strong experimental and theoretical in- 
terest to the half-metallic ferromagnetism in Heusler al- 
loys the number of the studies of the exchange interac- 
tions in Heusler alloys is still very small. The first contri- 
bution to the density functional theory of the exchange 
interactions in these systems was made in an early paper 
by Kubler et. al. , |25j where the microscopic mechanisms 
of the magnetism of Heusler alloys were discussed on the 
basis of the comparison of the ferromagnetic and antifer- 
romagnetic configurations of the Mn moments. Recently, 
the studies of the inter-atomic exchange interactions in 
several Heusler compounds were reported by the present 
authors and Kurtulus et al. 0000. 

The purpose of the present work is a first-principles 
study of the exchange interactions and the tempera- 
ture of magnetic phase transition for four semi Heusler 
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compounds NiCrZ (Z=P, Se, Te) and NiVAs that were 
recently predicted to be half-metals 0, 0- For all 
four systems we find the Curie temperature substantially 
exceeding room temperature. We discuss the relation 
between the value of the Curie temperature and half- 
metallicity. We demonstrate strong dependence of the 
effective exchange interaction between 3d atoms on the 
sp element (Z constituent). 

The paper is organized as follows. In Sec. II we present 
the calculational approach. Section III contains the re- 
sults of the calculations and discussion. Section IV gives 
the conclusions. 



TABLE I: Lattice parameters and magnetic moments (in /is) 
of NiVAs, NiCrZ (Z=P, Se, Te) and NiMnSb. 





a(A) 


Ni 


V,Cr,Mn 


Void 
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Cell 


NiVAs 


5.85 


-0.02 


2.05 


0.06 


-0.09 


2.0 






-0.03 


2.06 


0.06 


-0.09 


2.0 


NiCrP 


5.59 


-0.01 


3.08 


0.07 


-0.15 


3.0 




5-65(GaAs) 


-0.06 


3.16 


0.06 


-0.16 


3.0 


NiMnSb^ 


5.64 


0.19 


3.86 


0.04 


-0.10 


4.0 


NiCrSe 


5.64 


0.24 


3.64 


0.12 


-0.01 


4.0 




5-65(GaAs) 


0.23 


3.65 


0.12 


-0.01 


4.0 


NiCrTe 


5.84 


0.24 


3.68 


0.11 


-0.03 


4.0 




5-87(l n p) 


0.23 


3.70 


0.11 


-0.03 


4.0 



COMPUTATIONAL METHOD 
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The semi Heusler compounds crystallize in the C1&- 
type structure. The lattice consists of three interpen- 
etrating fee atomic sublattices. Compared with full 
Heusler alloys the forth atomic sublattice is vacant. In 
the calculations, the atomic positions of the fourth sub- 
lattice are occupied by empty spheres. 

The calculations are carried out with the augmented 
spherical waves method [3(J within the atomic-sphere 
approximation (ASA) |8ll | . The exchange-correlation po- 
tential is chosen in the generalized gradient approxi- 
mation j^. A dense Brillouin zone (BZ) sampling 
30 x 30 x 30 is used. 

For each compound we performed calculations for two 
values of the lattice parameter: the theoretical equilib- 
rium parameter [liJEl an d the lattice parameter of a bi- 
nary semiconductor ( Tabic GJ that can be considered as a 
possible substrate for growing the corresponding Heusler 
alloy: GaAs for NiCrP and NiCrSe and InP for NiVAs 
and NiCrTe. The radii of all atomic spheres are chosen 
equal. 

To calculate the interatomic exchange interactions we 
use the frozen-magnon technique |33l | and map the re- 
sults of the calculation of the total energy of the helical 
magnetic configurations 

s„ = (cos(qR n ) sin#, sin(qR„) sin#, cos#) (1) 

onto a classical Heisenberg Hamiltonian 

(2) 



eff 



where is an exchange interaction between two Cr (V) 
sites and Sj is the unit vector pointing in the direction of 
the magnetic moment at site i , R n are the lattice vectors, 
q is the wave vector of the helix, 9 polar angle giving 
the deviation of the moments from the z axis. Within 
the Heisenberg model the energy of frozen-magnon 
configurations can be represented in the form 



where Eq does not depend on q and J(q) is the the 
Fourier transform of the parameters of exchange inter- 
action between pairs of Cr(V) atoms: 



= £ J ° R ex P( ic l ' R )- 

R 



(4) 



Calculating E(9, q) for a regular q-mesh in the Bril- 
louin zone of the crystal and performing back Fourier 
transformation one gets exchange parameters Jor, be- 
tween pairs of Cr(V) atoms. 

The Curie temperature is estimated in the mean-field 
approximation (MFA) 
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and random phase approximation (RPA) 
1 1 
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M N ^ w(q) 



u/(q) (5) 



(6) 



E(9,q) = E (9)+ sin 2 9 J (q) 



(3) 



where w(q) = ^[J(0) — J(q)] is the energy of spin-wave 
excitations. 



RESULTS AND DISCUSSION 
DOS and magnetic moments 

The two lattice parameters used in the calculations re- 
sulted for all systems in very similar physical properties 
(see, e.g., Tables llTTHl . Therefore the most of the results 
will be presented for one lattice constant. All the discus- 
sion in this section is valid for both lattice spacings. 

In Fig. ^ we present the calculated electron densities 
of states (DOS) of the ferromagnetic states of the four 
Heusler compounds. All systems are found to be half- 
metallic with the Fermi level lying in the semiconducting 
gap of the minority-spin channel. Our DOS are in good 
agreement with the DOS presented in Refs.jU 0]. 
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FIG. 1: Spin-projected total density of states of NiVAs and 
NiCrZ (Z=P, Se, Te) 

In Table [I] the calculated magnetic moments are col- 
lected. Since the systems are half-metallic, the magnetic 
moments per formula unit are integer: for NiVAs, 

ZHB for NiCrP and for NiCrSe and NiCrTe. The ma- 
jor part of the magnetic moment comes from the second 
formula atom (V,Cr). Small induced magnetic moments 
are found on Ni and sp atoms. 

Exchange constants and Curie temperature 

Figure presents the frozen-magnon dispersion for 
one direction in the reciprocal space. Additionally to 
NiCrZ (Z=P, Se, Te) and NiVAs we present for compari- 
son the results for NiMnSb [2^| . The five systems can be 
subdivided into two groups. One group contains NiVAs, 
NiCrP and NiMnSb. Here the frozen-magnon dispersions 
are monotonous and resemble, visually, a simple cosinu- 
soid. The second group contains NiCrSe and NiCrTe 
and is characterized by nonmonotonous dispersions with 
a maximum close to the center of the q-interval (Fig. Et) ■ 
Note, that the sp elements (the third chemical-formula 
constituents) within each of the group belong to the same 
column of the Mendeleev's table whereas for different 
groups these columns are different. The importance of 
the valency of the sp element for magnetic properties of 
Heusler alloys has been already observed in our earlier 
studies 

The calculated exchange parameters are given in 
Fig. ^p. Since the inter-atomic exchange parameters are 
the Fourier transforms of the frozen-magnon dispersions 



they reflect the properties of the dispersions: The ex- 
change parameters belonging to the same group show 
similar qualitative behavior. On the other hand, there is 
strong difference between the systems belonging to dif- 
ferent groups (Fig. Et>) ■ In the first group (Fig. left 
panel), the strongest exchange interaction takes place be- 
tween nearest magnetic 3d atoms. This strongest interac- 
tion determines the cosinusoidal form of the correspond- 
ing magnon dispersion. The sizable interaction between 
the second- nearest magnetic 3d atoms describes the de- 
viation of the dispersion from a simple cosinusoid. 

In the second group of compounds (Fig. 03, right 
panel), the strongest interaction is between the second- 
nearest magnetic atoms. Because of the decreased period 
of the Fourier component corresponding to the second 
exchange parameter the dispersions are nonmonotonous 
and have the maximum not at the boundary of the Bril- 
louin zone but inside of it. 

A remarkable feature of the exchange interactions is 
their short range character: the leading contribution into 
the Curie temperature of all systems is provided by the 
interactions within the first two coordination spheres. 
The interactions with further coordination spheres are 
very weak and can be neglected in the calculation of the 
Curie temperature. The interaction between Ni atoms 
and the interaction of Ni with V and Cr are very weak 
and are not presented. 

The calculated exchange parameters are used to eval- 
uate the Curie temperature ( Table HT)l . It is important 
to note that the similarity of the form of the magnon 
dispersions within one group of compounds is not ac- 
companied by a quantitative closeness of the curves (Fig. 
0}). Therefore the Curie temperature can differ strongly 
for compounds belonging to the same group. 

The Curie temperatures are estimated within two dif- 
ferent schemes: MFA and RPA. For all systems and for 
both theoretical schemes the calculated Curie temper- 
atures are substantially higher than room temperature. 
The MFA always gives the value of T c that is larger than 
the corresponding RPA value and usually overestimates 
the experimental Curie temperature |37l l38| . The reason 
for the difference of the MFA and RPA estimations is 
a different weighting of the spin-wave excitations within 
two calculational approaches: The MFA takes all exci- 
tations with the same weight whereas the RPA gives a 
larger weight to the excitations with lower energy. The 
RPA weighting is better grounded from the viewpoint of 
statistical mechanics. 

The analysis of the calculational data allows to make a 
number of important conclusions. First, there is a strong 
influence of the sp atom on the value of T c . Indeed, the 
comparison of NiCrSe and NiCrTe that differ by the sp 
atom shows that the Curie temperature changes from 
about 500 K in NiCrSe to about 800 K in NiCrTe. This 
strong dependence reveals the sensitivity of the exchange 
interactions and the Curie temperature to the details of 
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FIG. 2: a)Frozen-magnon energies as a function of the wave vector q (in units of 2n/a) in NiVAs, NiMnSb and NiCrZ (Z=P, 
Se, Te) for V-V, Mn-Mn and Cr-Cr interactions, respectively. b)Interatomic exchange interactions in NiCrZ (Z=P, Se, Te) 
and NiVAs. The distances are given in the units of the lattice constant. The calculational data for NiVAs and NiCrZ are 
presented for the theoretical equilibrium lattice constant. The data for NiMnSb are from Ref. |2S|. 



the electron structure. 

An interesting feature of the calculated Curie temper- 
atures (Table ITTfl is a large difference between the MFA 
and RPA estimations for the first group of compounds 
in contrast to a small difference for the second group. 
Characterizing the relative difference of the MFA and 

RPA values of the Curie temperature by the relation 
{ T MFA _ T RPAy T RPA we get for the first group of com _ 

pounds a large value of 20-24% compared to a small value 
of 5-8% in the second group. This feature reflects the 
properties of the corresponding frozen-magnon spectra. 

In MFA, the Curie temperature is determined by 
an arithmetic average of the magnon energies while in 
RPA T c is determined by the harmonic average of the 
same quantities (the first average is always larger than 
the second). In terms of magnon energies, T^ IFA is 
equal to T^ PA in the case that the magnon spectrum 
is dispersion-less: the magnon energy does not depend 
on the wave vector q. 

From Fig. we see that the frozen-magnon curves of 
the second group of compounds are flat in the second half 
of the q interval demonstrating here a very weak disper- 
sion. On the other hand, the first group of compounds 
have considerable dispersion in this part of the q interval. 
Comparing now the low-q parts of the curves we notice 
that the curves of the second group of compounds lie 
higher than the first group. Therefore, the relative con- 
tribution of the low-energy magnons to the RPA value of 
the Curie temperature is smaller in the second group of 



compounds. 

This combination of features of the wave- vector depen- 
dencies of the frozen-magnon energies is responsible for 
a larger (smaller) difference between the RPA and MFA 
estimations of the Curie temperature of the first group 
(second group) of compounds in the two cases. 



Curie temperature and half-metallicity 

An important question concerning the magnetism of 
the half-metallic systems is the relation between half- 
metallicity and Curie temperature. Indeed, a number of 
studies has shown that the half-metallicity can stimulate 
an increase of the Curie temperature [sfl liH l4lj . The 
analysis of Fig. allows us to establish a correlation be- 
tween the value of the Curie temperature and the energy 
distance S between the Fermi level and the upper edge 
of the semiconducting gap. This quantity determines the 
spacing between the highest occupied spin-up state and 
the lowest empty spin-down state. For very small 5 of 
0.03 eV in NiCrSe we obtained the lowest Curie temper- 
ature of 508 K. On the other hand, for large S in NiCrP 
and NiMnSb we obtained the Curie temperature substan- 
tially above 800 K. 

Since the value of the Curie temperature is determined 
by the magnetic excitations to interpret the 5 — T c corre- 
lation we need to understand the origin of the influence 
of the S value on the spin-wave energies. 
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FIG. 3: Spin-projected total density of states of NiVAs, 
NiCrP, NiMnSb, NiCrSe and NiCrTe near Fermi level. Also 
shown are the RPA estimations of the Curie temperature and 
the energy distances between the Fermi level and the bottom 
of the spin-down conduction band. 



The magnon energies reflect the energy prize for the de- 
viations of the atomic moments from parallel directions 
pjil I35I [3(|. In the ground-state ferromagnetic config- 
uration the spin-projection is a good quantum number 
and the electron states with opposite spin projections do 
not interact. The deviation of the atomic moments from 
parallel directions leads to the mixing of the majority- 
spin and minority-spin states. The hybridization of a 
pair of the states leads to their repulsion. As a result, 
the energy of the lower state decreases (bonding state) 
and the energy of the upper state increases (antibonding 
state). If the lower state is occupied and the upper state 
is empty this process leads to the decrease of the energy 
of the magnon making ferromagnetic state energetically 
less favorable. The smaller is the energy distance between 
interacting states the stronger is the effect. 

Coming back to the half-metallic compounds consid- 
ered in the paper we note that the hybridizational in- 
teraction of this type takes place between the occupied 
majority-spin states below the Fermi level and the empty 
spin-down states at the bottom of the conduction band. 
The distance between these states is given by parameter 
8. Therefore this process provides a mechanism for the 
correlation between 8 and Curie temperature. 

Since the strength of the hybridizational repulsion in- 



TABLE II: MFA and RPA estimations of the Curie tempera- 
tures in NiVAs, NiCrZ (Z=P, Se, Te) and NiMnSb. The third 
column represent the relative difference of the MFA and RPA 
estimations. 
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lattice constant of GaAs 
lattice constant of InP 
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creases with decreasing energy distance between inter- 
acting states the negative contribution to the spin-wave 
energies is larger in the case of smaller 8. However, the 
correlation between parameter 8 and T c or, more general, 
between the half-metallicity and T c should not be consid- 
ered as a universal rule. The hybridizational repulsion 
considered above is only one of numerous processes aris- 
ing in a complex multi-band system with the deviation 
of the atomic moments from the parallel directions. The 
combined result of these processes cannot be predicted 
without the direct calculation of the excitation energies. 
Such a calculation must take into account the complex- 
ity of the electron structure of a real system. Indeed, 
the comparison of NiVAs and NiCrTe shows that both 
systems have the same 8 of 0.14eV. However, their Curie 
temperatures differ strongly. 

Also our numerical experiments with varying the lat- 
tice parameter in half-metallic systems have shown that 
in some systems the loss of the half-metallicity caused by 
this variation does not lead to the decrease of the Curie 
temperature |2sj . 



CONCLUSION 

We have systematically studied exchange interactions 
and Curie temperature in predicted half metallic semi 
Heusler compounds NiVAs and NiCrZ (Z= P, Se, Te). 
The calculations are performed within the parameter- 
free density functional theory. The RPA and MFA are 
used to estimate the Curie temperatures. We show that 
the behavior of exchange interactions in these systems 
vary strongly depending on the Z constituent. 

The exchange interactions are short range with the 
leading contribution to the Curie temperature provided 
by the interactions within the first two coordination 
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spheres of the magnetic 3d atoms. The predicted Curie 
temperatures of all four systems are considerably higher 
than room temperature. The relation between the half- 
metallicity and the value of the Curie temperature is dis- 
cussed. The combination of a high spin-polarization of 
charge carriers and a high Curie temperature makes these 
Heusler alloys interesting candidates for spintronics ap- 
plications. We hope that the present study will provide 
a guideline for experimental work stimulating the fabri- 
cation of these materials 
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